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Abstract 

The electronic response of a two-dimensional (2D) electron system represents 
a key quantity in discussing one-electron properties of electrons in semicon- 
ductor hetero junctions, on the surface of liquid helium and in copper-oxide 
planes of high-temperature superconductors. We here report an evaluation 
of the wave-vector and frequency dependent dynamical quadratic density- 
response function of a 2D electron gas (2DEG), within a self-consistent field 
approximation. We use this result to find the Z\ correction to the stopping 
power of a 2DEG for charged particles moving at a fixed distance from the 
plane of the 2D sheet, Z\ being the projectile charge. We reproduce, in the 
high-density limit, previous full nonlinear calculations of the stopping power 
of a 2DEG for slow antiprotons, and we go further to calculate the Z\ correc- 
tion to the stopping power of a 2DEG for a wide range of projectile velocities. 
Our results indicate that linear response calculations are, for all projectile 
velocities, less reliable in two dimensions than in three dimensions. 
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I. INTRODUCTION 

Since the pioneering work of Bohm and PinesS, the conduction electrons in a metal have 
been described as a three-dimensional (3D) gas of electrons in a neutralizing uniform positive 
charge. The dynamical linear density-response function of a 3D electron gas (3DEG) was 
evaluated by Lindhardi in the so-called random-phase approximation (RPA), in which each 
electron is assumed to move in the external field plus the induced field of all electrons. The 
wave-vector and frequency dependent dynamical quadratic density-response function of a 
3DEG has also been evaluated!, by going beyond linear response theory. The knowledge 
of this quantity has been proved to be of great importance in discussing the properties of 
electrons in a variety of 3D systems@, and, in particular, in explaining the experimentally 
observed difference between the electronic energy losses of protons and antiprotons moving 
through a solidii. 

The suggested existence of two-dimensional electron layers in metal-insulator- 
semiconductor structures and on the surface of liquid helium led several years ago to a 
great activity in the study of a two-dimensional electron gaj, where the electrons are con- 
fined to a plane and neutralized by an inert uniform rigid positive plane background. The 
2D electron system has also been considered in discussing the physics of new-class materials 
such as copper-oxide planes of high-temperature superconductors!. It has been found that 
electrons confined to a 2D layer show sometimes interesting properties not shared by a 3D 
electron system. For instance, for a 2D metal the plasma frequency goes to zero in the long- 
wavelength limit, in contrast to the 3D situation!. SternS evaluated the dynamical linear 
density-response function of a 2DEG in the RPA, and calculated the plasmon dispersion 
and the asymptotic screened Coulomb potential. Also, much effort has gone into studying 
the ground state energy0 and the excitation spectrum of a 2DEG§. The stopping power 
for a fast particle moving parallel to a 2DEG was first evaluated in the RPA, within linear 
response theory, by Horing et a£l, and the effect on this quantity of finite temperature0, 
local field corrections^ and recoilEl has also been considered. Nonlinear calculations of the 
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stopping power for slow protons and antiprotons have been performed only very recentlyEzil! 
on the basis of a scattering theory approach, the scattering cross sections being calculated 
for a statically screened potential. 

In this paper we present results for the dynamical quadratic electronic density-response 
of a 2DEG to longitudinal external fields of arbitrary wave-vector and frequency, which we 
evaluate on the same level of approximation as the RPA linear density-response function of 
StenM In order to illustrate the usefulness of the knowledge of this quantity, we consider, 
as an example, the stopping power of a 2DEG, and we use the quadratic density-response 
function to evaluate the Zf nonlinear correction to the stopping power of a 2DEG for particles 
of charge Z\e moving at a fixed distance from the 2D plasma. 



II. DYNAMICAL ELECTRONIC RESPONSE 

We consider a uniform electron system of density n at zero temperature. Linear and 
quadratic density-response functions of this system, x(x, x') and Y(x, x', x"), with x = (r, t), 
may be introduced in connection with the response to the presence of a time- varying external 
influence, after an expansion of the induced electron density in powers of the external poten- 
tial V(x). According to time-dependent perturbation theory, the induced electron density is 
given, up to second order in V(x), by (we use atomic units throughout, i.e., e 2 = h = m e = 1) 

n ind (x) = J dx' X (x,x')V(x') + J dx' J dx"Y(x,x' ,x")V(x')V(x"), (2.1) 

where 

X (x,x') = -ie(t-t') < * | \Ph(x),Ph(x , )\ |*o > (2.2) 

and 

Y(x,x',x") = &(t - t')Q(t f - t") < * | [[pH(x),p H (x')\rPH(x")\ 1*0 > /2 

+Q(t - t")G(t" - t') < tf | [[p H (x),p H (x")],p H (x')} |*o > /2. (2.3) 
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Here |^o > denotes the normalized ground state, pu is the density fluctuation operator 
Ph = Ph — n , where p H is the exact Heisenberg density operator in the unperturbed 
system, and Q(x) is the Heaviside step function accounting for causality. 

In a self-consistent field or random-phase approximation (RPA), it is assumed that the 
electron density induced by an external potential can be replaced by the electron density 
induced in a non- interacting electron gas by the sum of the external potential, V(x), and 
the potential created by the induced electron density itself, V md (x). Hence, in this approx- 
imation: 

n ind {x) = J dx' X °{x,x') [V(x') + V ind (x'j 

+ fdx' f dx"Y°(x, x', x") \v{x') + V md {x')] \v(x") + V ind (x")] . (2.4) 



Here x°( x , x> ) an d Y°(x, x', x") are 'free-particle' linear and quadratic density-response func- 
tions, and 

V md {x) = J dx'v{x,x')n ind (x'), (2.5) 

where v(x, x') represents the instantaneous Coulomb interaction. Introducing Eq. (2.5) into 
Eq. (2.4) and keeping in Eq. (2.4) only terms up to second order in V(x), one finds the 
following integral equations for RPA linear and quadratic density-response functions: 

X(*^-M*u«) + J^I^MWM4,*>) (2.6) 

and 

+J dx[ J dx'Wixi.x'Mx'^xVYixlx,,*,), (2.7) 

where K(x, x') is the so-called inverse dielectric function: 

K(x, x') = S(x -x') + J dx"v(x, x")x{x", x'). (2.8) 

The integral equations (2.6) and (2.7) have, within many-body perturbation theory, 
the simple diagrammatic interpretation shown in Fig. 1, where empty and full bubbles 



represent non-interacting and interacting linear density-response functions, x°( x ii x 2) and 
x(xi,x 2 ), respectively. Similarly, empty and full triangles represent non-interacting and 
interacting quadratic density-response functions, Y°(xi, x 2 , £3) and Y(x±, x 2 , £3), respec- 
tively, and dashed lines represent the Coulomb interaction, v(x,x f ). Thus, RPA linear and 
quadratic density-response functions are represented diagrammatically by summing over the 
infinite set of diagrams containing one (see Fig. la) and three (see Fig. lb) strings of empty 
bubbles, respectively. 

In the case of a homogeneous 2DEG, the electrons are free to move in two spatial di- 
mensions, having their motion constrained in the third dimension. Thus, assuming time 
invariance, we define the Fourier transforms 

Xq = J d 2 ri J dhe-^-^^-^xiruh-^h) (2.9) 

and 

Y qim = J d 2 ri y dhj d 2 r 2 J dt 2 e- i [ t n-( r i- r 2)-^( t i-*2)] e - i [( t i 1+t i 2 )'( r2 - r3 )-( a;i+a;2 )(* 2 - t3 )] 

xy(r 1 ,t 1 ;r 2 ,t 2 ;r 3 ,t3), (2.10) 

where ri, r 2 and r 3 represent two-dimensional position- vectors in the 2D plane, and q is the 
trimomentum q = (q, q°). Hence, within the RPA we find: 

X q = X° q + X° q v qXq (2.11) 

and 

= Yq,-qiKqiKq-qi + Xq V qYq,-qn (2-12) 

where 

K q = l + v qXq (2.13) 

and v q = 27r/|q|. 

For a non-interacting Fermi gas, the ground state is obtained by filling all the plane wave 
states inside the Fermi sphere of radius qF = \/2/r s , r s being the average interelectronic 



distance (rig 1 = 7rr 2 ). As in the case of a 3DEG0, we find non-interacting linear and 
quadratic density-response functions to be§ 



"X.q 



and 



Y 



o 



d 2 k 

(2^)2 



n k (l - n k+q ) 



1 - n k )^k+q 



(^k+q - Wk) + ir/ 



n k (l - rr k+q )(l - n k + qi ; 



d 2 k 

(2tt) 2 (g° + w k+q + ir/Xg? + u k - w k+qi + ir?) 
(1 - rr k )n k+q rr k+qi 

(g° + uj k - w k+q + ir])(q1 + ^ k - w k+qi + ir/) 
rr k+q (l - n k )(l - rr k+qi ) 



+ 



+ 



-g° + ^ k+q - a? k - ir/) (-g° + g? + u; k+q - cj k+qi - ir/) 

(1 ~ ^ k +q)^ k ^ k + qi 

-g° + u; k+q - a; k - ir/) (-g° + g? + ^ k+q - w k+qi - ir/) 

"k+gi(l - ^ k )(l - rr k+q ) 

-g? + ^k+ qi - Wk - ir/)(g° - g? + w k+qi - ^ k+q + ir/) 

(1 ~ rr k+qi )ra k ra k+q 

-gf + ^ k+qi - Wk - ir/)(g° - g? + w k+qi - ^ k+q + ir/) 



(2.14) 



+ 



?-?i)]> 



(2-15) 



where n q = 6(g^ — |q|), w k = k 2 /2, and r/ is a positive infinitesimal. 

Analytical evaluation of Eq. (2.12) results in the non-interacting linear density- response 
function of Stern0. As for the non-interacting quadratic density-response function, we first 
sum occupation numbers in Eq. (2.15) to find 



y: 







d 2 k 



1 



1 



r^k 



+ 



+ 



g u + ^ k - ^ k + q + ir/ gi + ujy, - w k+qi + ir/ 
1 1 



-g° + w k - w k+q - ir/ -(g° - g°) + w k - w k+q - qi - ir/ 
1 1 



-gf + ^ k - w k+qi - ir/ (g° - g?) + w k - u; k _( q _ qi ) + ir/ 



For the real part, we find§: 



(2.16) 



Re 



Y 



(Iq,qi + Iq,qi) (^-q,-Q+qi + I-q,-q+qi) + 



-31,9-91 1 -q\,q-q\i 



(2-17) 
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where 



1 



27r|q||qi| sin x 

Ai sin x 



{arctan 



ylsinx 



+ arctan 



A — A\ cos x 



Ai — A cos x 
— sgnAi arctan 



sgnA arctan 



A — A\ cos x 



sin X\l A 2 ~ Tf 
A\ — A cos x 



Q(A 2 -q 2 F ) 



(2.18) 



and 



I' 



9,91 



-e(q 2 F sm 2 x-G) 



1 



4|q| |qi| sin x 



(2.19) 



Here A = q°/\q\ - |q|/2, A l = g?/|qi| - |qi|/2, and G = ^A 2 - 2AAi cosx + ^i, X being 
the angle between q and qi. 

As for the imaginary part of Y q °-q 1 we define!!! the function H qm , which can be repre- 
sented in terms of a sum over hole and particle states and gives the second order contribution 
to the so-called absorption probability, as demonstrated in Ref. |6|. We find, 



Im 



y; 



9,-91 



H, 



9,91 



H, 



9i,9 



H, 



9-91, -91' 



(2.20) 



where 



§ 



H, 



9,91 



U9,9i 



/- 



9,-9+91 



+ 



q - Qi)] 



(2.21) 



and 



1 



27r|q||qi| sinx 



In 



Ai — A cos x + sin xy q\ — A 2 



Ai — A cos x — sin xy Tf ~ A 2 



(2.22) 



In particular, at low frequencies an expansion of H qm in powers of the frequency q° gives, 



after retaining only the first-order terms, 



H. 



9,91 



6(4g 



2 f-H 2 ) 



4(|q| cosx - |qi|) 



t q qi 



|q - qil 2 - 4^ sin x \/4g|-|q 



(2.23) 



In the static limit (q° — > 0) both first and second order contributions to the absorption 
probability, lmx q and H qqi , are proportional to the frequency q°, as in the case of a 3DEG. 
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III. ELECTRONIC STOPPING POWER 



We consider an ion of charge Z\ moving with constant velocity v at a fixed distance h 
from a 2DEG of density n . The Coulomb potential of this moving particle has the form 

V(r, z; t) =Z 1 \r-vt + {z- h)k\~\ (3.1) 

where r represents, as in Eqs. (2.9) and (2.10), a two-dimensional position- vector in the 
2D plane, and z denotes the coordinate normal to the 2DEG which we consider to be 
located at z = 0. Hence, in order to obtain the induced potential, we introduce this time- 
varying external potential into Eq. (2.1), and Eq. (2.1) into Eq. (2.5). We note that 
density-response functions of a 2DEG have their z arguments localized to the 2D plane by 
positional 5 functions, we Fourier transform, and find, up to second order in Z\. 

V ind (r,z;t) = Z 1 J ^e^-^HqKW+^x^ 

, 72 f d2t i f r i(q-r-^MIqlN + (lqil + lq-qiiy»l 7 . V ». ». (o O) 

J (27r) 2 J (27r) 2 t/ q 2 q,w;-qi,-c<; 1 t ; qit ; q-qi? 

where u = q ■ v and u)\ — qi ■ v. 

The stopping power of the 2DEG is simply the retarding force that the polarization 
charge distribution in the vicinity of the projectile exerts on the projectile itself^, and is 
given by 

S = — J d 2 r J &z5{y - wt)5{z - h)VV ind (r, z; t) ■ v. (3.3) 
Substituting Eq. (3.2) into Eq. (3.3), we havell 

"V / WT^ I p«- ,Mtl,llth - ,# '>K«,--l'VW. (3.4) 
In the RPA, \q an d Yq,-qi are obtained from Eqs. (2.11) and (2.12) or, equivalently, from 

VqXq = K q -l (3.5) 



S 



and 

Yq,-qi — K q Y®_ qi K qi K q _ Ql , (3.6) 

where 

K q =(l- X lv q y 1 . (3.7) 

The linear contribution to the stopping power of Eq. (3.4), which is proportional to 
Zf, was evaluated in the RPA by Horing et aJll at T = 0, and similar calculations were 
presented by Bret and DeutschB at finite temperature. On the other hand, the quadratic 
contribution to the stopping power of Eq. (3.4) is, in the RPA and for a geometry with the 
ion-beam in-plane within the 2D electron layer {h — 0), equivalent to the result derived in 
Ref. |5], within many-body perturbation theory, as the energy loss per unit path length of the 
projectile, the integration space being now in two dimensions instead of three dimensions. 
This contribution to the stopping power, which is proportional to Zf, discriminates between 
the energy loss of a proton and that of an antiproton, and appears as a consequence of losses 
to one- and two-step electronic excitations generated by both linearly and quadratically 
screened ion potentials, as discussed in Ref. ||. 

In the case of slow intruders (v — > 0), only the low-frequency form of the response 
enters in the evaluation of Eq. (3.4). At zero frequencies both linear and quadratic density- 
response functions are real, Im x q an d H q qi being at low frequencies proportional to the 
frequency q°. Thus, retaining only the first-order terms in the frequencies, both Z\ and 
Z\ contributions to the stopping power are found to be proportional to the velocity of the 
projectile, as in a 3DEG. For the RPA quadratic contribution to the stopping power we find, 
after insertion of Eq. (2.23) into Eq. (3.4), the following result: 



S 1 



(3) = 8vZf r . d|q| f°d| qi | r d xe -^\+M+\^ h {l±ll, (3.8) 



Jm-H |q-qil 



where 

fx = 0(2<7f - |q|)<o^ qi ,oiWo<o ;qi! o ( 3 - 9 ) 



and 



|q| (|q| cosx - 


Oil) 


7T 


Oil 


'|q- 


qi 


2 — Aq F sin 2 xj 



/2 =- y l 2 ^-|q|)— 7T, ; o ■ 9 \ ^q.O-Kqi.O-Kq-qi.O- (3.10) 



IV. RESULTS 

In the low- velocity limit — > 0), the stopping power can be evaluated to all orders in 
Zi, on the basis of a free-electron picture, with the additional assumption of independent, 
individual, elastic electron scattering, and it is easily found to be, in both two and three 
dimensions, proportional to the velocity of the projectileH By using density-functional 
theory (DFT)§ to calculate the self-consistent potential generated by a static charge sub- 
merged in a 3DEG, Echenique et evaluated the nonlinear stopping power of a 3DEG 
for slow ions. Nonlinear calculations for an in-plane projectile in 2D have been performed 
recently by using a nonlinearly screened scattering potential based on the Sjolander-Stott 
theoryEl, with the use of a nonlinear version of the linearized Fermi-Thomas potential^ in 
which the screening constant is determined from the Friedel sum rule, and on the basis of a 
self-consistent potential^ obtained within DFT as in Ref. |26| for the 3D 



case. 



The nonlinear calculations presented in Ref. [18] for the energy-loss of slow protons and 



antiprotons are very close to the calculations of Ref. \Z7\ when exchange and correlation (XC) 
contributions to the DFT scattering potential are excluded. These calculations are correct 
to all orders in Zi, and they represent, therefore, a good check for our quadratic response 
calculations, which should be exact in the high-density and/or low Z\ limits. Nevertheless, 



the approaches of Refs. [Ty|,[18] an< ^ E3 nave t ne limitation of being restricted to very low ion 
velocities [v << Vf, vf being the Fermi velocity), while our nonlinear Z\ corrections are 
valid for arbitrary non-relativistic velocities. 

The low-velocity limit of the quadratic stopping power, as obtained from Eq. (3.8) 
for an in-plane (h = 0) moving charge and divided by the velocity, is represented in Fig. 
2 by a solid line, as a function of the electron-density parameter r s . Stars and squares 
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represent the full nonlinear contribution to the stopping power for antiprotons reported in 



Refs. [18] and |27|, respectively, multiplied by a factor of —1, showing an excellent agreement, 
in the high-density limit, with our Zf nonlinear contribution. For r s > 2, higher order 
corrections become important, and for r s > 3, Z\ contributions (dashed line) are smaller 
than the Zf correction, indicating that the external potential cannot be treated, for these 
electron densities, as a small perturbation. In the case of a projectile moving at a given 
h distance above the 2D plasma the external perturbation is obviously diminished and, in 
particular, for h = l/qp the quadratic stopping power is smaller than the linear one for all 
electron densities (see the inset of Fig. 2). The full nonlinear contribution to the stopping 



power for protons reported in Refs. [Tl| and ^7] for h = is also represented in Fig. 2 by 
circles and crosses, respectively, showing large differences with our Zf calculations for all 
electron densities. These differences appear as a consequence of perturbation theory failing 
to describe electronic states bound to the proton, which in 2D systems can be supported by 
arbitrarily weak attractive potentials^. 

As the velocity of the projectile (v ) and/or the distance from the 2D plasma (h) increase, 
the ion potential becomes a relatively smaller perturbation and the Zf contribution to the 
stopping power for antiprotons may, therefore, be expected to approximately describe the 
full nonlinear contribution to the stopping power for arbitrary values of v and h, as long 
as r s < 2, and also for lower densities (r s > 2) as the velocity and the h distance increase. 
Substitution of the full RPA linear and quadratic 2DEG density-response functions of Eqs. 
(3.5) and (3.6) into Eq. (3.4) results in the quadratic stopping power plotted by a solid line 
in Fig. 3, as a function of the impinging projectile velocity, for r s = 1 and h = 0, and also 
for r s = 1 and h = 1/qF (see the inset of Fig. 3). The dotted line represents the low- velocity 
limit, as obtained from Eq. (3.8), and dashed lines represent linear RPA contributions to 
the stopping power of Eq. (3.4). 

The quadratic contribution to the stopping power of a 2DEG presents properties not 
shared by the 3DEG. First, the range of validity of the linear velocity dependence of this 
contribution to the stopping power (see Fig. 3) persists only up to velocities much smaller 
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than the Fermi velocity, in contrast to the 3D situation in which the linear velocity de- 
pendence persists up to velocities approaching the Fermi velocity!. Second, at velocities 
around the plasmon threshold velocity, for which the projectile has enough energy to excite 
a plasmon^, the ratio between Zf and Zf contributions to the stopping power increases, 
again in contrast to the 3D situation. Furthermore, we have found that the increase of 
this ratio at the plasmon threshold velocity becomes dramatic as the electron density de- 
creases. Of course, the ratio between quadratic and linear contributions to the stopping 
power decreases very rapidly with h (see the inset of Fig. 3), but the relative, and for large 
r s dramatic, increase of this ratio at the plasmon threshold velocity persists for all values 
of h that are not much larger than 1/qF- For a 2DEG the group velocity of the plasmon 
wave nearly coincides with the plasmon threshold velocity of the projectile (see Fig. 4); also, 
the inclusion of short-range correlations, which are ignored in the RPA, is known to have a 
substantial effect on the plasmon dispersion^. We interpret the anomalous enhancement of 
the ratio between quadratic and linear contributions to the stopping power at the plasmon 
threshold velocity and small electron densities (r s > 1) as a result of neglecting, within the 
RPA, short-range correlations between the electrons of a 2D system, since these correlations 
are non-negligible for all values of h as long as r s is not small. 



V. CONCLUSIONS 

We have presented an analytical evaluation of the wave-vector and frequency dependent 
non-interacting quadratic density-response function of a 2DEG. We have used this result to 
find, within a self-consistent field approximation, the Zf correction to the stopping power of 
a 2DEG for charged recoiless particles moving at a fixed distance from the 2D plasma sheet. 
We have reproduced, in the high-density limit, previous full nonlinear calculations of the 
stopping power of a 2DEG for slow antiprotons, and we have gone further to calculate the 
Zf correction to the stopping power for a wide range of projectile velocities. We have found 
that the Zf contribution to the stopping power of a 2DEG presents properties not shared 
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by the 3DEG. On the one hand, the range of validity of the linear velocity dependence of 
the Zf contribution to the stopping power persists only up to velocities much smaller than 
the Fermi velocity, and, on the other hand, an anomalous enhancement of the ratio between 
quadratic and linear contributions to the stopping power at the plasmon threshold velocity 
is found, within the RPA, at small electron densities (r s > 1). Also, our results indicate 
that linear response calculations are, for all projectile velocities, less reliable in 2D than in 
3D. 
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FIGURES 

FIG. 1. Diagrammatic interpretation of the RPA integral equations (2.6) and (2.7) for (a) linear 
and (b) quadratic density-response functions, respectively. Electron-hole empty and full bubbles 
represent non-interacting and interacting linear density-response functions, respectively. Empty 
and full triangles represent non-interacting and interacting quadratic density-response functions, 
respectively. The interacting RPA quadratic density-response function is obtained by summing 
over the infinite set of diagrams that combine three strings of empty bubbles (two-electron loops) 
through an empty three-electron loop. Dashed lines represent the electron-electron bare Coulomb 
interaction. 

FIG. 2. Low-velocity limit of the Z\ stopping power, as obtained from Eq. (3.8) (solid line) for 
h = and Z\ = 1, divided by the velocity, as a function of r s ; the corresponding Z\ stopping power 
is represented by a dashed line. Full nonlinear contributions to the stopping power for antiprotons 
\Z\ = — 1] (stars and squares) and protons \Z\ = 1] (circles ans crosses) have been obtained by 
subtracting RPA linear calculations from the results of Ref. [ll| (stars and circles) and Ref. |27] [with 
XC contributions to the DFT scattering potential excluded] (squares and crosses), and dividing by 
Z\\ thus, the negative values at r s > 1 simply mean that the full nonlinear stopping power lies, in 
the case of protons, below the linear results. The inset exhibits Z\ and Z\ contributions to the 
stopping power (dashed and solid lines, respectively) for h = 1/qF and Z\ = 1. 

FIG. 3. Full RPA Z\ stopping power, as obtained from Eq. (3.4) (solid line) for Z\ = 1, h = 
and r s = 1, as a function of the velocity of the projectile. The corresponding Z\ stopping power is 
represented by a dashed line, and the dotted line represents the low-velocity limit of the Z\ term. 
The inset shows the same results for h = 1/qF- 

FIG. 4. Two-dimensional RPA plasmon dispersion relation (solid line) and maximum energy 
transfer (co max = qv) at the plasmon threshold velocity (dotted line), as functions of the wave 
number. The electron density parameter has been taken to be r s = 1, thus the plasmon threshold 
velocity being v t h = 2.02i>o (t>o is Bohr' velocity). 
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Fig. 4 




